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The real time evolution and relaxation of expectation values of quantum fields and of quantum 
states are computed as initial value problems by implementing the dynamical renormalization group 
(DRG). Linear response is invoked to set up the renormalized initial value problem to study the 
dynamics of the expectation value of quantum fields. The perturbative solution of the equations of 
motion for the field expectation values of quantum fields as well as the evolution of quantum states 
features secular terms, namely terms that grow in time and invalidate the perturbative expansion 
for late times. The DRG provides a consistent framework to resum these secular terms and yields a 
uniform asymptotic expansion at long times. Several relevant cases are studied in detail, including 
those of threshold infrared divergences which appear in gauge theories at finite temperature and 
lead to anomalous relaxation. In these cases the DRG is shown to provide a resummation akin to 
Bloch-Nordsieck but directly in real time and that goes beyond the scope of Bloch-Nordsieck and 
Dyson resummations. The nature of the resummation program is discussed in several examples. 
D , The DRG provides a framework that is consistent, systematic and easy to implement to study the 

P-h ■ non-equilibrium relaxational dynamics directly in real time that does not rely on the concept of 

j. ' quasiparticle widths. 

PACS numbers: 11.10.-z,05.10.Cc,11.10.Wx,05.10.-a,02.60.Nm 
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I. 



INTRODUCTION 



The study of relaxation phenomena in quantum field theory or many body systems out of equilibrium is a timely 
topic of multidisciplinary interest with applications ranging from ultrarelativistic heavy ion collisions, to cosmology 
and condensed matter physics. The program of ultrarelativistic heavy ion collisions focuses on creating a new state 
of matter, the Quark-Gluon plasma[l|-|3| which existed about one microsecond after the Big Bang. At the ultrarel- 
ativistic heavy ion colliders RHIC at Brookhaven and the forthcoming LHC at CERN, this new state is conjectured 
to be formed and to exist for about 10 -22 seconds. Thus probing the physics of this short lived transient state 
experimentally requires the theoretical understanding of non-equilibrium phenomena on unprecedented short time 
scales. 

In cosmology the challenge of understanding non-equilibrium phenomena leading to thermalization and equilibra- 
tion, phase transitions and baryogenesisj^, ITfl ITU requires a novel set of tools to study relaxation in extreme 
environments and on very short time scales. 

In condensed matter, transport phenomena in mesoscopic svstems|l3j| an ultrafast spectroscopy in 
semiconductors [l4l 1 15j explore relaxation of non-equilibrium states on short time and distance scales. 

The usual tools to study relaxation and transport phenomena is based on kinetic equations which describe the 
dynamics of single (quasi) particle distribution functions or the concept of the damping rate of quasiparticles which 
is extracted from the imaginary part of the single particle self-energy and is associated with exponential decay. 

In gauge theories at high temperature it has been found that charged particles of hard momentum feature infrared 
singularities on their particle mass shells that lead to anomalous damping rates [l^. These infrared singularities are 
a consequence of the emission and absorption of transverse photons which are only dynamically screened by Landau 
damping. 

A study of the relaxation of charged excitations directly in real time by means of a Bloch-Nordsieck resummation 
of the propagator^^j revealed that this anomalous damping rate actually implies non-exponential relaxation. 

A novel approach, the dynamical renormalization group (DRG) was recently introduced |18l Il9j . It leads to a study 
of relaxation and kinetics directly in real time and reveals the non-exponential relaxation of charged particles in high 
temperature gauge theories. 

Th DRG approach to relaxation and kinetics borrows the main concepts from the study of the asymptotics in 
ordinary differential eauations poll2lll22]| . 

The dynamical renormalization group (DRG) has been recently implemented to study transport phenomena in high 
temperature ga uge theories: in particular the electrical conductivity has been extracted in high temperature quantum 
electrodynamics |23j . This approach to studying relaxation and transport phenomena directly in real time does not 
rely on the quasiparticle concept and instead extracts the real time dynamics from a resummation program that is 
akin to the renormalization group in critical phenomena and deep inelastic scattering. 

In this article we present the main concepts of the DRG explicitly obtaining the relaxation of expectation values in 
quantum field theory as well as the time evolution of quantum states. 

In section II and III we begin by setting up the fully renormalized initial value problem that describes the real 
time evolution of the expectation value of bosonic quantum field out equilibrium in linear response. In section IV we 
solve the real time evolution equations in several relevant cases by Laplace transform and highlight the case in which 
the emission and absorption of massless quanta introduce threshold infrared divergences. These divergences render 
the concept of a quasiparticle decay rate ambiguous and show that a Bloch-Nordsieck resummation of these infrared 
divergences lead to power law relaxation at zero temperature. At finite temperature the infrared divergences are more 
severe and a Bloch-Nordsieck resummation of the propagator is not readily available. 
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In section V we introduce the dynamical renormalization group resummation and compare to the results obtained 
by the Laplace transform solution highlighting the relevant resummations. This method allows to treat the most 
infrared severe cases of threshold divergences in a similar simple manner thus providing a very powerful resummation 
scheme. 

In section VI we present a scalar theory model that features the most relevant cases: stable particles and resonances 
and different cases of threshold infrared divergences which are common to gauge theories. 

In section VII the time evolution of states in quantum mechanics is solved in the DRG approach and a direct 
comparison with standard results on energy shifts and decay rates is established. 

In section VIII the resummations implied by the dynamical renormalization group, and the issue of separation of 
time scales are discussed in detail. Section IX summarizes the conclusions. 

Several appendices contain relevant mathematical details and a simple albeit illuminating example. 



II. FROM LINEAR RESPONSE TO AN INITIAL VALUE PROBLEM: 



We are interested in studying the real time evolution of expectation values of field operators. Consider a scalar field 
theory with an interacting Lagrangian density the expectation value of the scalar field <E> can be obtained from 
linear response to an external c-number source term J. The appropriate formulation of real time, non-equilibrium 
dynamics is that of Schwinger-Keldvsh |24| - |2fij in which a path integral along a contour in imaginary time is required 
to generate all of the non-equilibrium Green's functions. Along the contour the fields carry labels ± corresponding to 
the forward branch (+) along which the time variable takes values — oo < t < oo, and the backward branch (— ) along 
which t runs back from oo to — oo. The two branches describe the time evolution of a density matrix which requires 
the time evolution forward and backward in time (for details see refs.[3 HE HH)- 

The non-equilibrium Lagrangian density along this contour is therefore given bv|24j-|26[. 

C N eq[$ + , J] = £[$+] + J * + - £[*~] - J *~ . (ILL) 
The non-equilibrium expectation value of the scalar field in a linear response analysis is given by 

/+oo 

d 3 x' dt' G R (x - x', t - t') J{x', t') , (II.2) 
-oo 

with the retarded Green's function 

Gr(x -x',t- t') = [G>(x -x',t- t') -G < (x- x', t - t')] 6(t - = i(Mx, t), t'])e(t - t') , (II.3) 

where the expectation value is in the full interacting theory but with vanishing source. Consider an external source 
term that is adiabatically switched on in time from t — > — oo and of the form, 

J(f ,0 = J{x') e et ' Q(-t') , e -> 0+ . (II.4) 

The retarded nature of Gr(x — a? , t — t') results in that, 

tf>(g , t = Q) = 0o (f) , 4>{S, t < 0) = , (II.5) 

where (j)o(x) is determined by J{x) (or vice versa), and the second line follows in the limit e — > 0. For a given initial 
value (f>o(x) the current J{x)) can be found from from ea. (|II.2J) and the vanishing of the derivative for t < is a 
consequence of the retarded nature of Gr. 

The linear response problem with the initial conditions at t = given by eq. I|II.5|I can now be turned into an initial 
value problem for the equation of motion of the expectation value by using the (integro-) differential operator ^ 
inverse of Gr(x — of ,t — t') 

{s>t)( j)(x,t) = J(x,t) , 4>{S,t = Q)=4> Q {x) , 0(f,i<O) = O, (II.6) 



for the source term given by eq. (|II.4Jl . Within the non-equilibrium formulation, the equation of motion of the expec- 
tation value is obtained via the tadpole method[25l |26( and automatically leads to a retarded initial value problem 
by coupling an external source that satisfies eci. (|II.4|) . 
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III. A SCALAR THEORY 



Although the results obtained below are generic for bosonic, fermionic or gauge theories, we will focus our analysis 
on the interacting bosonic theory to highlight the main ideas in a definite setting. All of the steps can be generalized 
to more general theories, either fermionic, bosonic and or gauge theories. 

We now establish the general formulation for an arbitrary bosonic scalar field theory, deferring to specific examples 
in a later section. The type of theories under consideration are generally of the form 



£ = ^{d^f - - $ 2 + Jo $ + dm , (III.l) 



where the source coupled to the field <E> has been introduced to provide an initial value problem for its expectation 
value as explained in the previous section. 

The fully renormalizcd equations of motion are obtained by introducing the usual wave function, mass and coupling 
rcnormalizations, 

$ = $ r , J r = J , mlZ = m 2 r Z + Sm 2 , (HI.2) 

with a renormalization of the interaction terms included in the interaction Lagrangian. 

We now suppress the label r with all quantities being renormalized, and write the Lagrangian density in terms of 
renormalized quantities and counterterms 

£ = \{d^f - 1 TO 2 $ 2 + J $ + Cct + dnt 

dt = \ [{Z- 1) {d^f - m 2 $ 2 ] - ~ 6m 2 $ 2 . (III.3) 

The counterterms are adjusted in perturbation theory as usual, and the interaction Lagrangian is in terms of renor- 
malized quantities and the counterterms associated with the interaction. 
Writing, 

t) = ip(x, t) + <j>(x, t) with (ip(x, t)) = ; <f>(x, t) = ($(f , t)) , 
using the tadpole condition USES 

(V-(x,<))=0, (III.4) 

and taking spatial Fourier transforms, we find the following equation of motion linear in the amplitude for the 
expectation value for t > (when the external source J vanishes) 

/+oo 
ZZ(t-t')Mt')dt' = 0, (III.5) 
-oo 

with lo\ = k 2 + m 2 , and used that the adiabatic current vanishes for t > [see ea. (|II.4ll ]. Here — t') stands for 
the renormalized retarded self-energy. That is, — t') is obtained from the bare one after mass and wavefunction 
renormalization. 

As discussed in the previous section, the source is chosen so that <pk(t — 0) = <fik(0) 4>k{t < 0) = 0. £j?(f — t') is 
the retarded self-energy which is obtained systematically in perturbation theory. The Fourier representation of the 
retarded self-energy is given by, 

— e-'SiM, (III.6) 

-oo 27r 

It is very useful to express the self-energy in a spectral representation. In renormalizable field theories, usually one 
needs to make two subtractions in the dispersion relation. However, in the scalar models considered below, the wave 
function renormalization is finite and one subtraction is enough (we can set Z = 1). We can then write a spectral 
representation in the form, 

S fc ( w + iO) - S fc (0) = w / dw' , . Pk[UJ ' - . (111.7) 



— OO 



uj' (uj-uj' + iO) 
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where the spectral density Pfc(w) for real bosonic fields is an odd function of ui, 

Pfe(-w) = -pk{u) . (IHi 



Ea. (|III.7(l implies that pfc(w) is positive for positive uj. 
The real and imaginary parts of Efc(w) are given by 

/ + OC / \ r-\-oo / /\ 

dco' i P ^ L0) = 2u? V / du' - Pk ^ ] , Im£ fc M = -7T p fc (#J.9) 
-oo W ( w - W ) ./O W ( w - W ) 

where stands for the principal part. 

The retarded self-energy can thus be written in a spectral representation in the form 

/+oo r+oo 
Pk {uj) e - *"* du = -2 9(f) / duj Pk {uj) smujt . (111.10) 
-oo JO 



and the inverse Fourier transformation, 

f + OC 



E fc (w)=/ e iuJt Y,%(t) dt , lmw>0 
Jo 



(III.ll) 



In what follows we use the usual perturbative expansions in terms of a dimensionless coupling A for the self-energy 
and the renormalization counterterms 

E fc M - Asi 1) H + A 2 4 2) M + C(A 3 ) , Pfe H-A^ 1) H + A 2 pi 2) H + 0(A 3 ) . (111.12) 

The initial value problem with the condition that <fr(x, t < 0) = can be solved via Laplace transform. 

A. The Laplace transform 

Introducing the Laplace transform of 4>k(t) as 

/>oo 

(p k (s)= dt e~ st 4> k (t) , Res>0. (111.13) 
Jo 

and that of the retarded self-energy, 

/>oo 

E fe (s) = / dt e~ st Ef (t) , (III. 14) 

Jo 

which can be read off ea. ljIII.ll|) after the analytic continuation iui = — s, 

%k(s = -iw + 0) = E fc (w + iO) = Re E fc (w) + Um E fc (w) (III. 15) 

where the real and imaginary parts of the self-energy are given by eqs. H111.9|) ■ Notice that s approaching the imaginary 
axis from the right half-plane, s = —iu) + 0, corresponds to the causal choice u> + iO. 

The spectral representation of the Laplace transform E^(s) follows from ea . (|III.7|) by analytic continuation, 

E fc ( S )-E fc (0) = 2 S 2 (line) 



where we have used the property (|111.8fl in the last line and Efe(0) = Efe(0). 

Upon taking the Laplace transform, ea. (|III.5| ) becomes the following algebraic equation (see Appendix [XJ , 

[ S 2 +w 2 + E fc ( s )]^( s )- s fc (O) + ^^ [Efc(0) - Efc(s)] = 0. (111.17) 

The spectral density is consistently computed in a perturbative expansion in terms of the renormalized coupling 
constant. The renormalization aspects are well known and therefore we do not well into them. We just emphasize 
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that the initial value problem is fully renormalized formulated in terms of the fully renormalized self-energy Efc(ti> + zO) 
for which the relation (|III.15() holds in terms of the fully renormalized Sfe(s). 
The solution of Eq. (|111.17|l is therefore given by 



/ x Mo) 

Vk{s) = — — 



1 <4 + s fc (o) 

s 2 +uj 2 k + S fc (s) 



(111.18) 



This expression manifestly implies that the Laplace transform method to obtain the real time evolution is equivalent 
to a Dyson (geometric) resummation of self-energy insertions. 

The real time evolution is obtained from the inverse Laplace transform given by 

Mt) = J e st M*) ■ (111.19) 
where the contour T runs from — ioo to +ioo parallel to the imaginary axis in the Re s > half-plane. We note that 



s = is not a pole since the residue vanishes identically in ca. (|III.18JI . 

Notice that ea . 18|) provides a non-perturbative solution for the time evolution of the expectation value of the 
field for a given perturbative expression for Efe(s). Ea. ljIII.18() is an explicit non-perturbative resummation of the 
Schwinger-Dyson series. For a given order in perturbation theory for the self-energy this result corresponds to a 
geometric (Dyson) series with the self-energy calculated at a given order in perturbation theory. 

The information about the initial state is encoded in £fc(s). The self-energy is computed perturbatively with 
propagators containing the information about the initial state. For example, in the thermal case they depend on the 
initial temperature, and for non-zero density they depend on the chemical potential. 

In order to compute the integral ea. (|III.19|l it is convenient to deform the contour to the left s-plane. The singularities 
of Sfe(s) are usually along the imaginary s axis, such as poles at the particle dispersion relations and production and 
thermal cuts. 

IV. THE TIME EVOLUTION VIA LAPLACE TRANSFORM 

We now study in detail the solution of the time evolution via Laplace transform in several relevant cases. 

A. Isolated particle pole: stable particles 

We begin by studying the simple case in which the spectral density for to > has support above the position of the 
particle pole f2fc(~ u;^), namely 

Pk(uj) ^ for il k < uth < w < oo . (IV. 1) 

In this case s 2 + ui 2 + Sfe(s) vanishes at two points along the imaginary axis at s — ±iQk, with flk being the real 
solution of the equation 

nt - uj 2 k - z k (n k ) = o . (iv.2) 

Since pkfilk) — the imaginary part of the self-energy vanishes at this point. This situation corresponds to a stable 
particle. 

To lowest order in the coupling A we find, 

fife = LUk + SuJk , Su)k = tt~ ReEi^(wfe) . (IV.3) 

lu>k 

The Laplace transform (|III.19(1 is performed by wrapping the contour to the one shown in fig. ^where- the cuts along 
the imaginary axis run along u> t h < |Ims| < +oo with Lo t h being the particle production threshold. 

Deforming the contour T to the left yields the integral as the sum of the pole contributions plus the integral over 
the discontinuity of the cut along the imaginary s-axis. This cut may be a production cut for |Ims| > Uth or at 
finite temperature a thermal cut along \Ims\ < k (see section IVTj) below. In all these cases, provided the poles are 
separated from the cuts, which is the case for stable particles we can write the integral in ea. l|III.19|l as 



= ^ = ^ooB(n fc t) + [ w g + E fc (o)] ^ co f^ 

0fc(O) L J 7-oo w [c^-c^-ReSfcH] 2 + \bn£ k (w)Y 
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FIG. 1: Contour for the time evolution. 



where 



E fe (0) 



To lowest order in the coupling A we find, 

Z k = 1 - A 



1 ORcS fc (^) 
2fi fc 



w ' (uj 2 -uiD 2 



(IV.5) 



(IV.6) 



this wave function renormalization is finite irrespective of the choice of renormalization counterterm Z even in renor- 
malizable theories where yOfc(w) 



u) as u> 



Evaluating the expression ea. l|IV.4|) at t = we arrive at the sum rule 

r+ °° p*(o;) 



Z k = 1 



[«g + E*(0)] / °°- 

J-oo w 



ReE fc (w)r + [ImE fc (w) 



(IV.7) 



which is easily confirmed to lowest order in the coupling A by using the expression for Sfe given by ea. pil.9p . the 
expansions (|III.12|I and neglecting Sr^ in the denominator on the right hand side of eq. (|IV.7|) . 

The integral over the cut in ea. l|IV.4|) is usually too complicated to be computed in close form. However, we can 
evaluate its asymptotic behavior for late times provided we know the behavior of p(k, oj) near threshold. 

For example, in the zero temperature case the integral over the cut in ea. (|IV.4|l extends over \ui\ > ui t h where 
u)th stands for the first production threshold. Typical production cuts vanish at threshold with a power law at zero 
temperature, 



Pfc(w) = A (w - uith) 



(IV.8) 



where A is a positive constant. For example for the production cut of two massive bosonic particles, the threshold 
behavior is determined by a — i and uith = yj (2m) 2 + k 2 . 
We find from ea. (|IV.4| l for late times, 



S(k, i) (cut) 
xRe 



S fc (0)] 



dui 



Pfe(w) cos(wi) 



2 A [cj 2 



— e (w - ujth) 



w [ w a _ w 2 _ ReE^,^)] 2 + [ImE fc (w)] 2 

2A(/c 2 + m 2 )r(a + l) 1 
^ [k 2 +m 2 - u? th + RcZ k (Lj th )f t^ 1 



uJtht + — {a + 1) 



(IV.9) 



We see that the cut contribution decreases for late times with a power law. The smaller is the phase space near 
threshold (the larger is a) the faster decreases the cut contribution for late times|27j|. 

Notice that ea. (|IV.9|) holds provided the particle pole is separated from the threshold. That is for uith > ^k- 
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B. Resonances: 



In this case the imaginary part of the self energy does not vanish at uj k . The position of the (quasi) particle poles 
is determined by the equation 

s 2 p +wl + X k (s p )=0. (IV. 10) 

Writing s p = —i Q k + 7^ ; 7^ <C Q k with both Q, k and 7^ real, and comparing real and imaginary parts we find 

Pfe(^fe) 



Re£ fe (O fe ) = , 7 fc = -sign(7 fc ) -k 



2Qi 



(IV.H) 



The first equation is the same as in the previous case eq. (ITV~2)l . Since, 4r^ > unless there are instabilities, it is 
clear that there is no solution for the second equation. Thus there is no complex pole in the first (physical) Riemann 
sheet, corresponding to a resonance or quasi-particle. Since there are no poles in the first (physical) Riemann sheet 
in the complex s-plane, the Laplace transform is obtained by wrapping the contour around the imaginary axis with 
the result 



*(M)eM,K + ^(0)] 



dtv 



p k (uj) cos(wt) 



w [w 2 - uj 2 k - ReE fc (w)] 2 



(IV.12) 



In weak coupling, the integrand of (|IV.12(I features a sharp peak near lo 
eq. (IIV.12|) is of the Breit-Wi gner form, 



7T p k (L0) 



ReEfcHr + pmSfeH]' 



[ImS fe (w)] 
VL k , near this peak the integrand in 

(IV.13) 



with Z k given by ea. (|IV.5|l and 



7T Pk(^k) 



20* 



1 



1 dRcS k {u 



2i1 k Ouj 



(IV. 14) 



To lowest order in weak coupling, Q k is given by ea. l|IV.3|l and 

7T p k (u) k ) 
1 k = — h 

2 LU k 



0{\ 2 



(IV.15) 



Since these poles are in the unphysical Riemann sheet there is no pole contribution deforming the contour T in 
ea. piL19jl . However, for small coupling the poles in the second Riemann sheet are near the cut and the sharp peak 
of the resonance dominates the cut contribution. 



(IV.16) 



In such a case, from ea. (|III.19(l we find 

Mt) = Z k fc (O) e~ r » 1 cos(n k t) + A 2 ) , 



instead of ea. (|IV.4| l. Ea. (|IV.15|) exhibits the Lorentz contraction of the width (dilation of lifetime) with respect to 
that in the rest frame (k = 0). 



C. Threshold singularities and Bloch-Nordsieck resummation 



Let us consider the case in which the spectral density has a threshold precisely at lo = u k , vanishing linearly at 
threshold. A concrete realization corresponds to a particle of mass m ^ that emits or absorbs a massless particle 
(at zero temperature) with Lagrangian, 

L = - (<9$) 2 + - (dir) 2 - - M 2 $ 2 - I $ 2 7T (IV.17) 

which is a particular case of the specific examples analyzed in section (|VI|) . 
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The self-energy to one-loop takes here at zero temperature the form, 



(4tt) 2 



k 2 



logfl 



O 



= \m z 



k 2 



log 



A: 2 



OA 



where we identify the dimensionless coupling A 
ultraviolet logarithmic divergence. 



(2 7T? 



(IV.18) 

-, and have adjusted the mass counterterm to cancel an 



The wave function renormalization is UV finite but features a logarithmic infrared divergence at s 



-{k 2 



m 



That is. 



is divergent at s 2 



-{k 2 



This infrared divergence entails that there is no longer an isolated 



particle pole, which now becomes the beginning of a logarithmic branch cut as a consequence of the emission and 
absorption of massless quanta. 

The spectral density to lowest order in A is given by 



k 2 



Q(lo 2 



m 2 ) sign(w) 



lu 2 — k 2 

We find that the inverse propagator in the Laplace transform in ea. (|III.18|) is given by 



s 2 + k 2 + m 2 + E fe (s) = (s 2 + k 2 + m 2 ) 



1 + A 



k 2 



log 



OA 



(IV.19) 



(IV.20) 



The form of the inverse propagator given by ljiV.20JI indicates that in this case the logarithmic threshold behavior is 
associated with a wave function renormalization. 

This situation is common to a wide variety of theories, in particular in Quantum Electrodynamics where the ab- 
sorption and emission of massless photons introduces logarithmic threshold divergences in the electron propagator|28| 
which are enhanced at finite temperature|l7l Il9j|. The logarithmic infrared diver genc e signals the breakdown of the 
perturbative expansion near the particle mass shell. The Bloch-Nordsieck method 28] provides a resummation of the 
leading logarithmic contributions and is equivalent to a renormalization group improvement of the propagator and the 



self energy. It leads 28] to the following resumed form of the inverse propagator near the mass shell s + k 



■ m 



0, 



s z + k 2 + m 2 + S fc (s) -> (V + k 2 + m 2 ) 



2\l-\ 



(IV.21) 



The inverse Laplace transform with this resummed propagator [see ea. (|III.19| can be done by deforming the contour 
leading to the following integral representation, 



4>k(t) = tf> k (0) 



2 sin 7rA 



cos cot 



GS-0 



l-A 



(IV.22) 



While the integral can be computed exactly in terms of Bessel and Struve functions the asymptotic long time 
limit is dominated by the threshold behavior at the lower bound of integration. We find for mt 3> 1 , 



l + O 



1 

Wfe t 



(IV.23) 



Thus a Bloch-Nordsieck or renormalization group resummation of threshold infrared divergences leads to an asymptotic 
long time relaxation featuring an anomalous dimension. 

A perturbative expansion of eq. I|IV.23(I in powers of the coupling A leads to 



4>k(t) = 4> k (0) cos(w fe t) 



l-A lnwfci + A 2 , 



to k t 



(IV.24) 



the logarithmic secular term signals the breakdown of the perturbative expansion. The Bloch-Nordsieck resummation 
of the infrared threshold divergences in the propagator leads to relaxation in terms of a power law with an anomalous 
dimension A. We will see below that the logarithmic secular term in eq. PV.24| will emerge in the perturbative 
expansion in real time and that the dynamical renormalization group will provide a resummation of these secular 
terms leading to the same result as obtained from the Bloch-Nordsieck resummation, given by eq. ljIV.23() above. 

Anomalous power law relaxation has been found in condensed matter physics in models of mesosc~ansport[ia 
and in the X-ray absorption edge singularities in metals 29J where the spectral density also features logarithmic 
singularities at threshold [2!|, as well as in the relaxation of fermionic excitations in cold dense 
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1. Threshold divergences 



Consider the case in which the spectral density near to — tOk is of the form 



k 2 



+ b k 



9(w - w fc ) sign(w) , 



(IV.25) 



with bk ^ 0. This corresponds to the nonzero temperature case in the model (|IV.17(I . 

Following the steps outlined above, we now find that near the particle mass shell the inverse propagator is now 
given by, 



s z + k z + m z + E fc (a) = (s z + k z + m z ) 



1 + A 



,2 /„2 , ;„2 i m 2 

■In 



k 2 



Am bk In 



s 2 + k 2 + m 2 



0{\ 2 ) 



(IV.26) 

and we have chosen Z = 1. Whereas the logarithmic divergence near the particle mass shell (from the zero tem- 
perature contribution) can be interpreted in terms of a wave function renormalization and exponentiated via the 
Bloch-Nordsieck or renormalization group resummation, the nonzero temperature contribution bk leads to a stronger 
divergence, since not only the derivative of the inverse propagator but now the inverse propagator itself diverges 
at threshold. This strong divergence cannot be resummed either by the renormalization group or Bloch-Nordsieck 
resummation. 

As we shall see below, the dynamical renormalization group in real time provides in this case a resummation that 
extracts the asymptotic long time behavior. 



V. DYNAMICAL RENORMALIZATION GROUP: 

We obtain here the perturbative expansion of the field expectation value and its late time behavior exhibiting 
secular terms. Then, we resum the secular terms using the DRG. 

A. Perturbative Expansion and Secular Terms 



Expanding the Laplace transform ea. (|III.18(l in powers of A yields, 

s«( s )-e«(o) 



+ A 



s{s 2 +u; 2 ) {s 2 +w 2 ) 2 



0(\ 2 ) 



(V.l) 



The real time evolution follows by inverse Laplace transform order by order in perturbation theory. It is convenient 
first to express the self-energy as a dispersion relation [see ea . 161 ] yielding, 



¥>*(*) = MO) 



S 2 + IjJ 2 



A 



s4 1} (0) 



2.1 r 

Jo w 



G(A 2 ) 



(S 2 +LU 2 ) 2 k J LU (S 2 + OJ 2 ) 2 (S 2 + LU 2 ) 

The inverse Laplace transform of ea. l|V.2|) is easily obtained comparing with the elementary integrals, 

dt e~ st cosw/jt , 

t dt e~ st sinukt , 



(V.2) 



s 2 + J 1 



k JO 
S 1 

S 



1 



( S 2 +u; z ) 2 ( S 2 +u 1 2 ) (u; 2 -u; z ) 2 J Q 
We find comparing eas. (|V.2| ') and (|V.3() and then using ea. ijIII.16(l . 



dt e st [cos uit — cos uikt ■ 



2u)k 



4>k(t) = (f> k (0) I cosu! k t + A 



dw p ( k 1} (uj) 



t 

2 u k 



t sinwfet] 



2wt V I — . r r. K v ""^,» (cos cut - cosLu k t) - — — sinuj k t ReE^MWfe) 

'o UJ (UJ 2 - LJ%) 2 > ' 



(V.3) 



0(X 2 ) 



(V.4) 
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The second term of order A turns to be the frequency shift SuJk given by eq. (jiV.3(l 
To first order in A we finally have, 



<f>k(t) = <fik{0) \ cos LUkt - t 6uk sinwfct + 2\Lu k V — 



diu p[} ] (w) 



(coswi — coswfci) + C(A 2 ) 



(V.5) 



The second term is a secular term: it grows as t. Furthermore, the frequency integral in the third term may also 
feature secular terms at long times. Its t dependence for late times is determined by the behavior of [ui) for 
u> — > LUk as we analyze below in several relevant cases. 

B. Dynamical renormalization group 

In perturbation theory we find that the solution of the evolution equation has in general the asymptotic form, 

+ c.c. (V.6) 



1 + A S { h 1} (t) + A 2 5< 2) (t) + non secular + C(A 3 ) 



f 1 2) 

with ' (t) containing only secular terms, namely terms that grow in time. 

The dynamical renormalization group absorbs the secular terms into a renormalization of the complex amplitude 
at an arbitrary time scale t. This is achieved by introducing a multiplicative renormalization in the following manner, 



<M0) = Mr) Rk(r) 
with i?fc(r) found systematically in an expansion in the coupling A, namely 

i? fc (r) = l + Ar[ 1) (r) + A 2 rf(r) + 0(X 3 ) 



(V.7) 



(V.8) 



The coefficients rl^ (r), which play the role of counterterms, are chosen to cancel the secular terms at the time t = r 
order by order in the perturbative expansion. To lowest order in A we find, 



The solution renormalized in this way at a scale r is to lowest order in A given by 



Mt)=Mr){l + \[si}\t)-Si 1 \r) 



non secular 



G(A 2 )} 



c.c. 



Since r is an arbitrary scale, the solution must not depend on it. We therefore request 

d 



(V.9) 
(V.10) 

Cv.li) 



This is the dynamical renormalization group equation. 

Inserting ea. HV.10|) into ea. ijV.llfl to lowest order in the coupling, the dynamical renormalization group equation 
becomes, 



ln[^(r)]=5«(r) 



with solution 



(pk(r) = (f>k(ro) e k 



(V.12) 



(V.13) 



Since 4>k{t) is independent of r thanks to the DRG ea. HV.ll|) . we can choose the arbitrary scale r to coincide with t. 
This yields, 



^ : i A' 2 . 1 



(V.14) 



This is the long time asymptotic behavior given by dynamical renormalization group. The overall multiplicative 
normalization can be adjusted up to O (A 2 ) using ea. (|V.6(l . The next sections extract the leading secular term Sj^'(t) 
in several cases. 
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C. Case I: Particles and quasiparticles 

We consider first the case where the frequency uj k is away from the multiparticle thresholds, namely 

Wfe ^ uj t h ; Pk(u > 0) ^ for uj t h < u < +oo . (V.15) 

That is, the spectral density is continuous, finite and non-zero at uj = u>k- This case treats on equal footing the 
situations corresponding to a stable particle uj k < ioth and that of a quasiparticle uj k > <-^>th- 

The long time limit t — > oo can be extracted from ea. l|V.5|) by using formulae IjB.lfl in appendix We find the 
asymptotic long time behavior to be given in this case by 



+ 0{\ 2 ) } + c.c. . (V.16) 



The perturbative solution displays secular terms that in this case grow linearly in time. The secular terms invalidate 
the perturbative expansion for late times t 3> — r ■ 

Clearly, a resummation scheme must be invoked to resum the perturbative expansion. 

In the notation of ea. HV.6(l we have, 

4°(t) = ^r \iTteE£\u, k )-wp£\u, k j\ . (V.17) 



k 



Therefore, ea. l|V.14|) for the DRG solution in the long time limit takes here the form, 



This can also be written as, 



2lul. 2u)i, 



"*(') = 5,, <>,(()) coMS),,./)e- rfct + 0(i A 2 ) • (V.I!)) 



with fife, Zfc and given by equations (|1V.3|) . (|1V.6|I and i|lV.15[l respectively. 

Comparing this DRG solution with that obtained from the Laplace transform, given by ea. l|IV.16|) . we see that in 
this case the dynamical renormalization group provides a resummation akin to that of the Dyson series. 

If pfe(wfc) = 0, namely the (stable) particle pole is below the multiparticle threshold, then — and the solution 
describes the asymptotic time evolution of a stable particle. On the other hand if u) k is above the multiparticle 
threshold, then the pole moves off the physical sheet into the second Riemann sheet becoming a quasiparticle pole. 
This case describes a decaying particle or resonance. 

D. Threshold singularities: 

We consider in this subsection the case where the spectral density is discontinuous at uj = ui% (case II) or it exhibits 
a logarithmic singularity at this point (case III). 

1. Case II: a discontinuous density at uj = ui k . 

Let us consider the case in which the spectral density around u> = uj k is given by equation (|IV.25|I 

This case reflects a discontinuity at threshold since p\ 
threshold uj ~ uj k the spectral density is approximated by 



This case reflects a discontinuity at threshold since p^ k \uj = uj k — 0) = O;/?^ (cj = u>k + 0) = A k b k and near 



p { k \uj) = A k [uj-uj k + b k ] [1 + 0(u - io k )} 0(w - w fc ) , A k ,b k >0. (V.20) 
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We can now extract the asymptotic long time behavior for <f>{t) given by equation (|V.5(1 by using the results given by 
eqs. (|B~21) . and |{B~T|) in appendix B. We find, 



fc ( f ) *^£» e to»f <j i + < , 



5w fc - A ^ fc &fc log(w fc £ e 7 x ) 



2uji 



tb k 



2b k 



logK t e 7 " 1 ) 



0(t°)\ +c.c + 0(X 2 ) , (V.21) 



i we see that the real linear secular term is given by r<pt with r 



With the convention 0(0) 

In particular, for bk = we identify the secular term —A ln.[u>fc£] which is manifest in the perturbative expansion 
ea. (|IV.24|) of the Bloch-Nordsieck resummed result given by ea. (|IV.23|) . 
Up to 0(A) we now find ea. lV.21|l given by ea. l|V.6|) with the secular term, 



Si 1) (t)=it 



Sto k - A Ak bk \og{uj k t e 7 x ) 



2 LUk 



tb k +[l 



2bk 



log(w fc t e 7 x ) 



The DRG resummation ea. ljV.14fl takes here the form, 



(V.22) 



4>k(t) *=°°0*(O) cos[S fc («)t] e- Tkt 




3k(t), T fe and z/fc are given by 



s*(t) 



w fc + 5ujk 



2uj k 



1 



XA k b k 

2 LLl k 

2bk 



ln^fcie 7 - 1 ) , 
irXA k b k 



r fe = 



4wt 



(V.23) 



(V.24) 



Thus, we see that the logarithmic mass shell singularity proportional to bk in the perturbative expansion of the 
propagator, ea. (|iV.26(l . is manifest in the real-time perturbative expansion given by ea. (|V.22|) as a secular term of 
the form tln[t]. The dynamical renormalization group resums this secular term into a logarithmic phase in the field 
expectation value ea. (|V.23|) . 

Notice from eas. l|V.20|l and ljIV.19|) t hat Ak = 2u k and bk = for the model (|IV. 1 7|> . Hen ce, we see in ea. (|V.24j) 
that Vk — X in this case. Therefore, eqs. JELll and lfv"23lt fully agree. In addition, ea. l)IV.23|) implies to = 2e 1 ju>k 
for the reference time scale. 

Thus, an important conclusion resulting from the comparison between the time evolution obtained from the Bloch- 
Nordsieck resummed propagator given by eq. ljiV.23(l for bk — and the dynamical renormalization group resummation 
in this case, is that the DRG resummation is a Bloch-Nordsieck resummation in real time. 

Moreover, while a Bloch-Nordsieck resummation of the Euclidean propagator for bk ^ given by eq. ljiV.26(l is not 
available, the dynamical renormalization group is able to resum the case b k ^ 0. 



2. Case III: a logarithmically divergent spectral density at lu — LUk- 

We now consider the case in which the spectral density is logarithmically divergent at LUk but does not feature a 
normal threshold at this value of frequency, namely 



C k 



log 



U) — U>k 



0(uj - LU k ) 



c k >o 



(V.25) 



In this case, using formulae l|B.3|l - (|B.4|) from appendix |BJ we find that the field expectation value takes the form of 
ea. (|V.6(l where the secular term Sj^'(t) is given by 



t 8wk — 



Xn C k 



In (/j t 



A 7T Ck 
2 LO k 



t ln^ie 7 " 1 ] 



(V.26) 
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Inserting ea. l|V.26|) in ea. (|V.14|l we find the DRG asymptotic long time behavior, 




. fc (0) cos[fi*(t)] e 



— 7fe t In [/it] 



+ C.C. , 



(V.27) 



with 



(lO k +Su! k )t 



ln(pt) 



(V.28) 



and 



Ik = 



irXC k 
2 uj k 



M = A* e7 1 > <^fc = o — ReS^(wfc) . 



(V.29) 



E. Summary of dynamical renormalization group resummation to lowest order 



We are now in position to summarize the main results obtained to lowest order in the dynamical renormalization 
group resummation program. 

The nature of the secular terms is completely determined by the behavior of the spectral density p k {Lo) near lo = uj k . 
The different cases studied above to lowest order in the coupling A describe fairly general behavior of the spectral 
densities in quantum field theories at zero and non-zero temperature. 

The important aspects of the results of this study can be highlighted as follows 

• i): when the spectral density p k (u>) vanishes at lo = co k the particle is stable and the only effect of the interaction 
is to renormalize the mass and wave function. This corresponds to the spectral density of the particle featuring 
an isolated pole away from the branch cuts associated with the multiparticle production. The resummation 
implied by the dynamical renormalization group is equivalent to the Dyson resummation of the propagator that 
leads to the renormalized dispersion relation for the particle. When p k (lo = Lu k ) ^ but Lo k > LOth with ujth being 
the multiparticle thresholds, namely the beginning and end points of branch cuts in the spectral density of the 
particle, the particle is unstable and its propagator features a pole in an unphysical sheet in the complex lo plane, 
namely a resonance. In perturbation theory this resonance is very close to the real axis and the propagator is 
of the Breit-Wigner form. The dynamical renormalization group in this case provides a resummation which is 
equivalent to the Dyson series. 

• ii): When the position of the multiparticle thresholds coincide with to kl infrared divergences arise near the 
particle mass-shell. This corresponds to the emission and absorption of massless particles, a situation which is 
common in gauge theories and in theories with spontaneous symmetry breaking with Goldstone bosons. The 
secular terms that determine the relaxation of the mean field in perturbation theory depend on the behavior of the 
spectral density near thresholds. The dynamical renormalization group resummation in the case p k (io k ± 0) = 
is similar to the Bloch-Nordsieck resummation of infrared divergences but directly in real time. No Bloch- 
Nordsieck resummation is available for the case II, which we solve by DRG and where the spectral density is 
discontinuous (p k (L0 k + 0) ^ 0). Furthermore, the DRG provides a consistent resumation in the case III with a 
logarithmically divergent spectral density at lo — to k . 



The secular terms associated with the different behaviors of the spectral density for these generic cases are sum- 
marized in the table below. 



Case 


Phi") = 


Secular term (t) 


DRG 

resummed amplitude 
<5(fc, t) 


Rcsummation 

done 
by the DRG 


I 


p k (uj = Uk) ^ ; p k {uj = uj k ± 0) ^ 
and continuous, ujk ^ <^th 






it(u, k +6u, k )-t ' P 2 fc i" fc) 

fc- fc 


Dyson 


II 


A k [uJ - LOk + b k ]<d(uj - u>k) 
A k , b k >0 


^ [ReS { k 1] (w fc ) - A fe 6 fe In (w fc t e^ 1 ) 

2"l \ w fe y V K ) 4iu k 


life) 


x e 2 ^fc 


Dyson + 

Bloch-Nordsieck 
for the kernel 


III 


C fc ln " 
C k > 




it(w fc +5a)fe) — 5-§^ t log(pt) 

i 7rC fc 

X (/2t) a k 


Dyson + 

Bloch-Nordsieck 
for the kernel 
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VI. SPECIFIC EXAMPLES IN A SCALAR THEORY: 



To highlight the main concepts in a simpler setting we study a scalar theory of three interacting fields $^ ; i = 1,2, 3. 
The Lagrangian density is given by 



3 r 



i 



-{d^f - - ml, d>? 



- 9o $1 $2 $3 + Jo $1 



(VI.1) 



where the source coupled to the field $1 has been introduced to provide an initial value problem for its expectation 
value as explained in section II, of course we could study simultaneously the initial value problem for the expectation 
values of all three fields, but our focus is to highlight the most important aspects in a simpler setting. 

The renormalization is performed by introducing wave functions, mass and coupling renormalizations as follows 



$ 4 = Zl' 2 $ i>r ; 



7 _ 7 i/2 T 



m j Zi = m r i Zi + 5m i 



; g Q y/Zi Z 2 Z 3 = g r Z g , i = 1,2,3. (VI.2) 



We now suppress the label r with all quantities being renormalized, and write the Lagrangian density in terms of 
renormalized quantities and counterterms, 



£ = E 

i=l 

c ct = ± 



1 



1 



mi 



1 



- (Zi - 1) (d^f - $2 _ Sm 2 $? 



-g(Z g - 1) $1 $ 2 $3 



(VI.3) 



The counterterms are adjusted in perturbation theory as usual. 
Some cases which are relevant for the discussion that follows are: 

• i: with mi > m% + m^, in this case the decay and recombination processes 1 <-> 2 + 3 are allowed. This will 
lead to the presence of a resonance, namely the quasiparticle pole will be in the continuum resulting in a decay 
rate. 

• ii: 771 1 = 7712 7^ ; 7773 = 0, this case will feature anomalous logarithmic relaxation as a consequence of the 
absorption and emission of soft massless quanta. 

In what follows, we will treat both cases by setting 7773 = and 7771 = 777 to present the main results in a simpler 
and more clear manner. 

The calculation of the self-energy of particle 1 at finite temperature is best performed in the imaginary time 
formulation. The Matsubara Green's functions for the fields are written as spectral representations in the form 



Gi(k,u n 
with the spectral densities 
Pi(po,k) 



1 



k.'i 



dpo ; — ; uj n 

PO - lU n 



2nnT 



1,2,3. 



2wfc 



[6(p - Uk,i) - 8(po +Wfc,i)] , Lul a = k 2 + m 2 , 7 = 1,2,3 



The one- loop self-energy for the particle 1 is shown in figure [21 and is given by 
£i,fc(iu;„) = ~g 2 T ^ 



d 3 p , , P2(P2,P) P3(P3,W+k\ 
dp2 dp 3 



(27T) 



P2 ~ W m p 3 - ILOr, 



(VI.4) 



(VI.5) 



(VI.6) 



The sum over the Matsubara frequencies can be done using the methods described in [^J leading to the following 
spectral representation for the self-energy 



Si,fc(io; n ) - Si k(0) = iuj„ 



Lu'(iuj n — uj') 



(VI.7) 
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FIG. 2: One loop contribution to the self-energy of particle 1. 



with the spectral density pk{uj) given by 
2 / d 3 P 



(2nf 



dp 2 dp 3 p 2 (p2,p) P3(P3, \P+k\) [n(p 2 ) - n(p 3 )} S(p 3 -p 2 -uj) ; n{p t ) = — 



1 



e t — 1 



i= 1,2,3 

(vi.s; 



Comparing this dispersive representation for the self-energy with that of equation 1111.7(1 we find that 

E fc (w) = Xi, fc ^„ = w + iO) (VI.9) 
A straightforward but lengthy calculation with m 2 = M ^ 0; m 3 — leads to the following spectral density 



Pk{oj) = A m 2 



w 2 - Wf T 



1 - e"" 1 - e~ v 
l-e u ~ 1 - e"~ 



1 - e~ u 1-e 
1 - e~ u ~ 1 - e" 



9(fc 2 - w 2 ) } sign(cj) 



e (w 2 _ w ^ + 



with 



2T(HTfe) 



M 
T 



W fe 2 = fc 2 + M 2 



(VI. 10) 



(VI.11) 



We identify the dimensionless coupling constant A = ( 2 nm)' 1 • 

The contribution to the spectral density for w 2 > I'K 2 corresponds to the two particle cut of one massive and one 
massless particle, and remains in the zero temperature limit. The contribution for uj 2 < k 2 corresponds to the Landau 
damping cut, only present at finite temperature and below the light cone. 

This spectral density models all of the situations studied in the previous sections: 

• Case I: to 7^ M . In this case u>k is away from the two particle threshold at Wk and above the light cone, thus 
away from the Landau damping region. If to < M the particle is stable and the self-energy correction amounts 
to a renormalization of the dispersion relation. If to > M now uik lies in the two particle continuum, and the 
particle pole moves off the physical sheet, in this case the particle becomes a quasiparticle with a finite width 
oc A. In both cases pk(y> — ^k) ^ and pk{^) ^ for uj ^ ujk- This situation describes the case I analyzed in 
detail above and corresponds to the most common cases of stable particles or resonances. 

• Case II: to = M. In this case the two particle threshold is precisely at uJk, but the two particle and Landau 
damping thresholds are well separated for any finite momentum fc. The behavior of the spectral density near 
the threshold at u>k is given by 



Pfe(w) = A TO 



— {uj 

TO 



UJk) 



I 1 " 



UJk 


+ k 




UJk 


- k 





6(w - ujk)- 



(VI. 12) 



This spectral density is of the same form as that of case II given by ea. HV.20|) . This case also applies to self 
energies of charged particles in gauge theories where the emission or absorption of massless gauge particles 
imply that ujf. is actually the beginning of a branch cut in the spectral representation of the charged particle 
propagator. While the (instantaneous) Coulomb interaction is screened at finite temperature by the polarization 
of the medium, transverse gluons (in perturbation theory) and photons are only dynamically screened by Landau 
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Thus the self-energy of charged particles in general gauge theories features the threshold 
infrared divergences associated with case II and described by the threshold behavior given by eas. l|V.20(l and 

• Case III: k >> m = M, or M = m = 0. In the case M — m and for a hard particle k >> m or alternatively 
for M = m = 0, the behavior of the spectral density ea. ljVI.10fl above for uik — > k is seen to be of the form 

(VI. 13) 

where fi(k) — 2T ^1 — e^^J. In this case the Landau damping cut merges with the two particle cut since 

LUk — k and the spectral density features a finite temperature logarithmic divergence at the position of the 
particle mass shell. This is the situation described by case III analyzed above and arises in the case of a hard 
charged particle, either fermionic or bosonic in high temperature gauge theories, where again the emission and 
absorption of transverse gauge fields, which are only dynamically screened introduces logarithmic divergences. 
In this case the notion of a damping rate which is associated with exponential relaxation is not the correct one 
to describe the relaxation of single quasiparticle excitations. 



- Am — In 
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VII. TIME EVOLUTION IN QUANTUM MECHANICS 



While we have focused the discussion above on the evolution of the expectation value of quantum fields, the DRG 
approach also offers a systematic method to study the evolution of quantum states. In particular, as explained below 
in detail, the DRG offers an alternative to the adiabatic switching-on procedure to generate exact eigenstates of the 
interacting Hamiltonian from the eigenstates of the non-interacting system. Along the way this DRG formulation 
extracts the energy shifts, wave function renormalization constants and offers a clear description of Fermi's Golden 
rule. 

Consider the time evolution of a quantum state in an interacting theory with Hamiltonian 

H = H + \H! , (VII.l) 

namely 

|V>(t)) = e - iff ('-*o) tyfo)) = e-^ U(t, t ) e m ^m )) . (VII.2) 
Where we have introduced the unitary time evolution operator in the interaction picture 

t t t' 

U(t,t ) =l—i\ I H!(t')dt'~X 2 j I Hi(t') Hi(t") dt' dt" + £>(A 3 ) , U(t ,t ) = 1 , (VII.3) 

with the interaction Hamiltonian in the interaction picture 

H^t) = e lHat Hi e- lHot . (VILA) 
It is convenient to pass to the interaction picture in which the states are given by 

m))i = e iff °* \m) , (vii.5) 

and their time evolution is given by 

m))i = U(t,t )\iP(t ))i. (VII.6) 

Clearly, the time evolution of quantum states is an initial value problem. 

The adiabatic, or Gell-Mann-Low theorem asserts that introducing the adiabatic time evolution operator U e (t,to) 
by replacing 

H z (t) -» e- £ l*l H!(t) , (VII.7) 
in the time evolution operator U(t,to) given by (|V11.3fl . the states 

= U £ (0, -oo) |n) , (VII.8) 
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constructed out of the eigenstates |n) of the non-interacting Hamiltonian Hq with energy E n , are exact eigenstates 
of the total interacting Hamiltonian H , and the exact eigenvalues are then given bv|33| 



E n + AE n — 



QPn|g|ttn) 
(*»!*») 



(VII.9) 



We now show how the dynamical renormalization group provides an illuminating alternative to this procedure. To 
make this approach more clear and to establish contact with the familiar results, we now consider the case in which 
\ip(t(j))i is an eigenstate of the non-interacting Hamiltonian Ho, namely 

|V(io))i = \n) ; H \n) = E n \n) . (VII.10) 

To highlight the main ideas of the DRG, we will focus on the evolution of the (persistence) amplitude, 

c n {t) = Hv(t))« , (mil) 

however the time evolution of the off-diagonal overlaps C m (t) — (m\i/)(t))i; m ^ n can be studied along the same lines. 
Inserting the identity J2 m \ m )( m \ = 1 appropriately, we find for the amplitude C n (t) the following expression, 



C n {t) = {n\U{t,t Q )i>{t Q )) i = C n {t ) 

'to Jt m ^ n 



X 2 

1 — i X (t — to) {n\H!\n) - — (t - to) 2 (n|i7/|n) 2 



-A 



Jtn Jtn l- 



(VII. 12) 



where we have written C n (to) — (n\ip(to))i despite our choice of initial state for which C n (to) = 1, to emphasize that 
the initial amplitude factors out. 

It proves convenient to introduce the spectral density 



Pn (w)= ]T \(n\Hj\m)\ 2 S(cu - E rl 



(VII. 13) 



in terms of which and after the straightforward integration over the time variables t',t", ea. i|VII.12|l becomes 



C n (t) = C n (t )< 1 -i A (t -t ) (n|£fi|n) - — {t - t Q ) 2 (^H^n) 2 



(VII. 14) 



A 2 / duo p n (ui) 



E n - UJ 



t-to 



sm[(E n - uj)(t - t )}\ l-cos[(E n -u)(t-t )} 



E n - LU 



o(x 3 



Using formulae IB. Ill and (|B.5|) of appendix |Bj we find the asymptotic long time limit to be given by, 

C n (t)= C n (t )\l-iX(t-t ) (n|ff z |n>-^ (t - t a ) 2 (n\H j\n) 2 - i(t - t ) X 2 ^ ^ U ' " 



E n — E n 



-A 2 7T (t - t ) Pn(E ri ) - X 2 ^ 



(E n — E m ) 2 \t — to 

where ^ m refers to the sum over all the states with E m ^ E n . We write the expression above as, 



(VII. 15) 



C n (t) — C n (to) 



1 + A S« (t) + X 2 S n 2 ^ (t) X 2 + O (j^-, A 3 ) 



(VII. 16) 



where X^S^(t) are secular terms (linear and quadratic in time respectively), the term —X Z n ' in eq. (|VII.16jl is the 
time independent term in ea. l|VII.15|) . 

Applying here the dynamical renormalization group, we introduce the amplitude renormalization in the form 



C«(to) = C„(r) Rn(r) ; R n (r) = 1 + A r^(r) + A 2 r n 2 \r) + 0(X 3 ) 



(VII. 17) 
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Hence, 

C n (t)=C n (T){l + \[rW(r)+SW(t) 



+ A 2 



f )(r) +S®(t) + rf'(r) S«(t)] - A 2 4 2) + <?(A 3 )} (VII. 18) 



The counterterms ?n (r) are required to cancel the secular terms at a time scale t = t, leading to 

r«(r) = -5W(r) , r«(r) = -5j?>(r) + [r«(r)] 2 . 



(VII. 19) 



The independence of the amplitude C n (t) on the arbitrary time scale r, namely dC n (t)/dT = leads to the dynamical 
renormalization group equation 



a 



X(r){ 



1 + A 



W(r)+5W(t)] +0(X 3 )}+C n (T){Xr^(r)+X" [f< 2 > (r) + (r) 5« (t)] + 0(A 3 )} =0 

(VII.20) 

where the dot stands for derivative with respect to r. Keeping terms up to second order in A yields, 



Cu{t) 
C»(t) 



?i 2) (r)-i (r r W(r)) : 



0(A 3 ) 



Notice that the ^-dependent pieces cancelled out identically. Such cancellation is necessary for the consistency of the 
DRG to second order in A . We find using the counterterms given by ea. l|VII.19|) . 



C n (r) = C n (t ) e x s nH T )+^ 2 {^(r)^ 



(VII.21) 



Using now the explicit form of the secular terms which are read off from eq. (|VII.15|) , we finally find the DRG solution 
to second order to be given by 



C n (r) = C„(Q) e 



with the energy shift and width up to 0(A 2 ) given by 



(VH.22) 



(VII.23) 



Notice that the terms in r 2 cancel identically in the exponent of ea. l|VII.21JI and l|VII.22(l . Again, such cancellation 
is necessary for the consistency of the DRG. 

Finally, inserting ea. ilVII.22fl into ea. (|VII.18|] and choosing the arbitrary renormalization scale r to coincide with 
the time t, we find the asymptotic long time behavior of the dynamical renormalization group resummed amplitude 
as 



C n (t) ^ Z n e-^ e-^ 



(VH.24) 



The energy shift AE n is clearly the same obtained in familiar perturbation theory, the decay rate r„ coincides with 
the transition probability per unit time obtained from Fermi's Golden rule, and the wave function renormalization 



i-a 2 E' 



(n\Hi\m)\ 2 d 



[E n — E m ) 



dE r . 



(E n + AE n 



(VII.25) 



is the same as obtained in usual perturbation theory and determines the overlap between the unperturbed state and 
the exact eigenstate of the Hamiltonian 33] . Thus, we see that the dynamical renormalization group leads to an 
alternative formulation of the construction of the exact eigenstates which allows to extract the energy shifts, widths 
and weights and that coincides with the familiar setting of quantum mechanics. 

There is an important byproduct of this exercise: by considering the perturbative expansion up to second order in 
the interaction Hamiltonian, we found secular terms that contain the square of the first order secular terms. These 
must cancel consistently and systematically in the final DRG equation since they are accounted for in the resummation 
furnished by the DRG. Indeed we found such cancellation in eas. (|VII.21|) and (I VII. 2 21) . This is the equivalent of the 
renormalizability of the perturbative expansion. Thus, this example not only provides a pedagogical framework to 
explore and confirm the dynamical renormalization group, but also manifestly shows the renormalizability in the sense 
that higher order terms that appear in the perturbative expansion, which are associated with the expansion of first 
order terms in the solution, cancel systematically in the dynamical renormalization group equation. 
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VIII. DYNAMICAL RENORMALIZATION GROUP: RESUMMATIONS 



The examples studied in detail above clearly indicate that the dynamical renormalization group provides a resum- 
mation of the perturbative series. Thus the following question arises : which type of resummation is implemented by 
the DRG?. The answer to this important question can be gleaned from the examples corresponding to the different 
cases discussed in the previous sections. 

• Case I: in the case in which the position of the particle mass shell near Uk is far away from the multiparticle 
thresholds in the perturbative expansion, then the DRG provides a resummation which is akin to the geometric 
Dyson sum of the propagator with the perturbative self-energy but keeping only the secular terms at each order. 
The DRG yields the asymptotic long time behavior of the expectation value of the field in the case of stable 
particles or unstable quasiparticles with a lifetime (width or decay rate) determined by the spectral density at 
the position of the mass shell of the particle. 

• Cases II and III: In the case in which there are threshold divergences in the propagator corresponding to the 
case when the position of the particle mass shell at u> ~ u>k coincides with the beginning (or end point) of a branch 
cut, the DRG provides a resummation of the self-energy kernel along with a resummation of the perturbative 
series for the propagator a la Dyson. That this is indeed the case can be understood from the study in section 
IIV Cl where Bloch-Nordsieck resummation of the propagator exponentiates the logarithmic threshold divergences. 
This resummation of the propagator is akin to that via the usual Euclidean renormalization group and leads 
to a power law relaxation, just as the scaling behavior of renormalization group improved correlation functions. 
The case of threshold divergences at finite temperature is much more severe by the enhancement of the infrared 
behavior by the Bose-Einstein occupation factors. In this case there is no readily available Bloch-Nordsieck 
or renormalization group resummation for the Euclidean propagator, but the DRG furnishes the resummation 
directly in real time. Thus, in conclusion, in all of these cases in which there are infrared divergences, typically 
arising from the emission and absorption of massless quanta, the DRG provides a resummation of the self-energy 
kernel along with the Dyson type resummation of the perturbative series. 

It is important to highlight that the DRG does all of these resummations automatically. Whereas in the usual 
approach, there is first the resummation of the Dyson series of the propagator, and in the case of threshold 
divergences the self-energy has to be further resummed independently, either by the Bloch-Nordsieck method 
or alternatively by the Euclidean renormalization group. Thus, insofar as extracting the asymptotic long time 
behavior of the expectation values, the DRG provides a powerful, systematic method with a straightforward 
implementation which furnishes several stages of resummation within a single approach. 



Secular terms and separation of time scales: For any finite time the perturbative expansion is well defined 
and is free of any divergences (ultraviolet divergences had been properly renormalized) . It is for long times that 
secular terms that grow in time invalidate the perturbative expansion. Consider the example of case I where the 
solution to first order in the coupling is given by ea. l|V.16l) . It is clear from this expression that the perturbation 
expansion breaks down at a time scale tb when 



A t b 



iReS^K) 



7T p\.'{Uk) 



2u k 



2uJk 



1 



(VIII. 1) 



This time scale is obviously of the same order as the relaxation time tb — t re i = T k where Tk is the width given 
by ea. HV.18(l . Secular terms dominate the integral in eq.{V3J) for sufficiently late times. Let t p be such time scale, 



namely, for t > t p we can approximate the integral in ea. (|V.5(l by the corresponding secular terms. As is clear from 
the explicit examples discussed above, the nature of the secular terms as well as the scale t p depend on the behavior 
of the spectral density near uj — ujk- The detailed behavior of the spectral density depends both on the field theory 
studied as well as the properties of the medium or the initial state chosen. 

Thus, the perturbative expansion features secular terms but it is still reliable in a time interval such that 



t„ <t <C t re i 



(VIII. 2) 



Therefore, if there is a wide separation of time scales such that t re i 3> t p , then, there is a wide region in time during 
which secular terms can be clearly identified and yet the perturbative expansion is valid. The secular terms can be 
extracted and renormalized during this large interval and the DRG resummation extends the result uniformly for 
t 3> t re i , namely beyond the validity of the original perturbative expansion. 

Since the spectral densities do not depend on the coupling (to lowest order) but the relaxation time scale is such 
that t re i ~ l/[m A] a wide separation of time scales is justified in the weak coupling limit. In higher orders the width 
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is always at least one more power of the coupling as compared to the spectral density. Hence, if there is a separation 
of scales in lowest order, such separation will remain in higher orders. 

Thus, the conclusion of this discussion is that the DRG allows a systematic improvement of the perturbative series, 
leading to a uniform expansion, in a weakly coupled theory in which a wide separation between the typical time scale 
beyond which there emerge secular terms and the relaxation time scale is justified. 



DRG and Euclidean RG: 



The connection with the more familiar Euclidean renormalization group can be established with the following 
example. Consider a theory of a single scalar field <p with a quartic self-interaction £j„ t = — ^y- 4 in four space-time 
dimensions defined with an upper momentum cutoff A. The 2^2 particle scattering amplitude with all equal 
external momentum is given by 



T^(p,p,p,p) = \ -^ X 2 t , t = In(-). (VIII.3) 



Thus we see that the usual logarithmic divergences in the scattering amplitude can be interpreted as secular terms in 
the perturbative expansion. 

Introducing Ah[k], the coupling renormalized at a scale k as 

A = X R [ K ] Z x [k] ; Z x [k} = 1 + X rZi [k] + 0(X 2 ) , (VIII.4) 

and choosing Z\[k] = hi(^) = r to cancel the secular term in the variable t in the scattering amplitude at the scale 
p = k, we find to second order in the renormalized coupling 



r (4) (p,P,P,p) = Xr[k] + | X 2 r [k] (t - 1) , i = ln(-) , r = ln(£). (VIII.5) 



3 ,2 r..i ^ , M\ _ ,_ / A \ 

■p 

Since the scattering amplitude cannot depend on the arbitrary renormalization scale «, namely 

d 



dT r^( P ,p,p,p) = o, (viii.6) 

the renormalization group equation for the coupling follows, namely 

Xr + Ix r = O^0 x = -^ X 2 r ; p x = k ^X r [k] . (VIII.7) 

Solving this renormalization group equation, where [3\ is recognized as the usual renormalization group beta function, 
and choosing the scale r to coincide with t, we obtain the renormalization group improved scattering amplitude 

r< 4 ) (p, P , P , P ) = Xr[k] . (viii.8) 

1 + | X r [k] ln(*) 

This expression reveals that the renormalization group is resumming the geometric series of one loop bubbles but only 
keeping the secular terms in the loops, namely the logarithmic terms. This resummation is tantamount to summing 
the leading logarithms. 



IX. SUMMARY AND CONCLUSIONS 



In this article we studied the real time evolution and relaxation of expectation values of quantum fields as well as 
the evolution of the states in quantum mechanics by implementing the dynamical renormalization group. 

The time evolution of expectation values of a bosonic quantum field is studied as an initial value problem. The fully 
renormalized equation of motion for the expectation value is obtained in linear response to an external current that 
is adiabatically switched on. We first studied the time evolution in terms of the Laplace transform of the equation 
of motion with a one-loop retarded self-energy, which is tantamount to a Dyson resummation (geometric series) of 
the self-energy. The real time evolution was obtained in several cases of interest: stable particles and quasiparticles 
(resonances) and in the case in which the propagator features threshold infrared divergences, which are a consequence 
of absorption and emission of massless quanta at zero or non-zero temperature. These threshold divergences invalidate 
the perturbative expansion and are ubiquitous in gauge theories. A Bloch-Nordsieck resummation leads to anomalous 
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relaxation in the case of logarithmic infrared singularities but is not readily available in the case of more severe infrared 
singularities. 

A strict perturbative expansion of the solution of the equations of motion l|III.5t features secular terms, namely 
terms that grow in time and invalidate the perturbative expansion. We straightforwardly obtain the perturbative 
solution of the equations of motion by expanding the Laplace transform solution cas. lIII.18j l- lIII.19jl in powers of A. 
One can alternatively solve by iteration the evolution equations (jIII.5j) Tsj], 

We then implemented the dynamical renormalization group (DRG) to resum these secular terms in all cases. The 
DRG effectively resums the self-energy kernel of the evolution equations (jIII.5jl while the Laplace transform 
eas. ljIII.18jl - ljIII.19p provides the exact solution with the input kernel (which is to first order in A here.) The DRG 
method provides a consistent and systematic approach that at once resums the secular terms. The DRG solution 
ca. (jV.14l) is a uniform asymptotic expansion for late times. 

The secular terms are linear in time in the case of stable particles or quasiparticles and logarithmic in time or linear 
times logarithmic in the case of threshold infrared divergences. In the case of stable particles and quasiparticles (reso- 
nances) the resummation implied by the DRG is the same as that of the Dyson series or Laplace transform. However, 
in the case of threshold infrared divergences the DRG provides at once a Bloch-Nordsieck type of resummation of the 
self-energy kernel and a Dyson series resummation of the solution. 

We then studied the time evolution in quantum mechanics as an initial value problem for the states by imple- 
menting the dynamical renormalization group. Again, the perturbative expansion features secular terms, which upon 
resummation via the DRG allow to extract the energy shift and width of quantum levels as well as the wave function 
renormalizations. Thus, the DRG provides a framework to systematically obtain energy shifts and width of quantum 
states. The consistency of the DRG to the second order of perturbations is proved here. Such consistency is analogous 
to the renormalizability of quantum field theory for the renormalization group. 

While we discussed primarily the case of bosonic theories, the method can be simply generalized to fermionic and 
gauge theories ^ lilt 

There are many distinct advantages to this method that make it a very powerful way to study relaxation directly 
in real time. While the most often used description of relaxation is based on the concept of quasiparticle widths 
and is associated with exponential relaxation, the DRG does not bias the description: if the secular terms are linear 
in time, indeed the concept of a quasiparticle width is the relevant one. On the other hand if the secular terms 
are more complicated, for example logarithmic or linear in time times logarithmic, the concept of a quasiparticle 
width and exponential relaxation is not the proper one. The dynamical renormalization group leads to the correct 
relaxational dynamics either power laws in the case of logarithmic secular terms or more complicated exponentials in 
the case of linear times logarithmic. These relaxational dynamics simply cannot be extracted from the usual concept 
of quasiparticle widths which are infrared diver gent 1 1 6l Il7| . 

Thus, one of the main conclusions of this article is that the dynamical renormalization group provides a powerful, 
systematic and easy to implement framework to study the relaxational dynamics of expectation values and the time 
evolution of quantum states directly in real time. This method implements a resummation of the perturbative 
expansion which yields a uniform asymptotic long time behavior and reveals the correct relaxational dynamics, all 
within one simple method. 

While we have illustrated the main aspects of the DRG to lowest order in the perturbative expansion, the next step 
is to extend the DRG calculations to higher loops in field theory. 
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APPENDIX A: RENORMALIZED EQUATION OF MOTION AND ITS LAPLACE TRANSFORM 

The wave function and mass renormalization [ea. (jIII.3jl ] are defined by the counterterms 5m 2 and Z which are 
consistently computed in a perturbative expansion in terms of the coupling constant. Hence, the renormalized 
retarded self-energy is related to the bare one by 

E fc (fl) = (Z - l)(s 2 + wjj) + S Mare (s) + Sm 2 , (A.l) 

Using this renormalized self-energy the tadpole condition ea. ljIII.4j) leads to the equation of motion (jIII.5jl . 
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We perform now the Laplace transform of ea. (|III.5| . We find, 



dt 



4> k (t) = s 2 <p k (*)-fa(P)-*4>k(P) 



(A.2) 



In the Laplace transform of the kernel term we change the integration variable A' by x = t — t' > and then the 
integration variable t by y = t — x. We obtain, 



dt e~ st 



dt' E fe (t - t') cj) k {t') = die 
Jo 



dx E fe (x) 4> k (t - x) 



e~ sx V k (x) dy e- sy 4> k (y) = dx e~ sx E k (x) 

J-x Jo 

E fc (0)-E fc (s) + Efe(s) (fik(s) , 



</> fc (o)+/ dy e~ s y My) 



(A.3) 



where we used that 4>k(t) = 0&(O) for t < 0. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF SPECTRAL INTEGRALS 

We summarize in this appendix the late time behavior of integrals over the density of states used in Sections III. 
In the formulas below p(y) stands for a smooth function for < y < oo. p(0) and p'(0) are finite. For large y ,p(y) 
are assumed to fall so that the integrals over y converge at infinity. 
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(B.l) 
(B.2) 
(B.3) 
(B.4) 
(B.5) 

(B.6) 



o 2/ 



sin yA 



dy 



p(y) = i p(0) [ln(/x A) + 7 - 1] + A / - [p(») - p(0) - y)] + - ) (B.7) 

y j Jo y \ r 



where A is a fixed positive number and 7 = 0.5772157 ... is Euler's constant and V stands for the principal part. 

Notice that the formulas are independent of the scale /i as one can easily see since the derivative with respect to 
fi of the r. h. s. identically vanishes. The scale /i has been introduced just to have a dimensionless argument in the 
logs. 

APPENDIX C: THE DYNAMICAL RENORMALIZATION GROUP IN A SOLVABLE EXAMPLE 

In this appendix we present the dynamical renormalization group through a simple example chosen for pedagogical 
reasons which allows to illustrate the fundamental features within a simple setting. 
We consider a damped harmonic oscillator: 



y + y = -ey , e < 1 



(C.l) 
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This equation admits as exact solution, 



y(t) = a e ' 1 cos ( t \j 1 - — + b 



(C.2) 



where the constants a and b are determined by the initial conditions. 
We seek a solution as a perturbative expansion in e, of the form 

y(t) = y (t) + e Vl (t) + e 2 y 2 (t) + 0(e 3 ) 

where the yt(t) are solutions to the following hierarchy of equations: 

Vo (t) + Vo (t) = , y\ (t) + yi (t) = -y (t) , y 2 it) + y 2 (t) = - yi (t) 

These equations can be solved iteratively by starting from the zero order solution 

y (t) = A e lt + c.c , 



using the retarded Green's function 



which is the solution of differential equation 



G ret (i) = B(t) sint 



dt 2 



G ret (t)=5(t) 



Up to second order in e, the perturbative expansion of the solution is given by 



y(t) = Ae l 



e e 2 2 e 2 ' 

1 — t-\ — r + i — t 
2 8 8 



c.c. + non-secular 



(C.3) 



Note that this solution contains secular terms that grow in t, the terms denoted by non-secular remain finite at all 
times. We see that this asymptotic expansion is only valid for et <g; 1. 

The dynamical renormalization is achieved by introducing the renormalization constant Z(t) in the form 



with the expansion 



A = A(t) Z(t) 



Z{t) = 1 + e Zi(t) + e 2 z 2 (t) + C(e 3 ) 



The coefficients Zi(r) are chosen to cancel the secular terms at a the time scale t = r. This is similar to choosing the 
wave function renormalization in field theory to absorb the UV (or IR) divergences at a given momentum scale. Up 
to C(e 2 ) we find 



*X(T) - ~ 2 



T T 

Zo(t) = i — 

K ' 8 8 



After this renormalization the solution 1|C.4|I takes the form, 

:2 



y{t,r) = A{t) e 



1 - - (t - t) + - (i - t) z + i - (t - t) 



c.c. + nonsecular 



(C.4) 



We have now an expansion in powers of e\t — r|. This asymptotic expansion is expected to be valid for e\t — t\ <C 1. 
That is, renormalization has improved the validity region of the solution from the neighborhood of t = [see ea. (|C.4ll ] 
to the neighborhood of t = t which is an arbitrary point. 

Since r is an arbitrary scale, the solution must not depend on it. We therefore request 



dy(t,r) 



= 



(C.5) 
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This condition leads to the following equation to this order 

2£U(!-«S)«„-.. 

where we have neglected terms of order higher than e 2 . Ea. ljC.6| is the dynamical renormalization group equation to 
this order. Ea. ljC.6|l can be easily solved with the solution 



A(t) = A(0) e -i T+l ^ T . 
Since the solution y(t, r) is independent of r thanks to ea. (|C.5fl we finally set t = t which yields, 

y(t) = A(0) e-t'+^-T)* + c.c. (C.7) 

Comparing eq. (|C.7|) with eq. (|C.2fl shows that eq. (|C.7(1 provides an uniform approximation to the exact solution l|C.2|) 
up to 0(e 3 ). That is, the dynamical renormalization group solution (|C.7|I differs from the exact solution (|C.2(1 by 
C(e 3 ) whatever is t. 

The solution l|C.7|) is nothing but an envelope of the family of solutions i|C.4() parametrized by t[23 • In other words, 
the r-family (|C.4() provides a solution for t near r while ea. (|C.7|) is a good approximation valid for all t. 
Further examples can be found in refs. 
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